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Motivated by recent transport and magnetization measurements in BSSCO samples 
[B. Khaykovich et at, Phys. Rev B 61, R9261 (2000)], we present a simple macroscopic model 
describing effects of inhomogeneous current distribution and shear in a layered superconductor. Pa- 
rameters of the model are deduced from a microscopic calculation. Our model accounts for the strong 
current non-linearities and the re-entrant temperature dependence observed in the experiment. 



I. INTRODUCTION 



Transport measurements are widely used in studies of 
vortex dynamics of high-T c superconductors. When the 
current distribution in the sample is not homogeneous, 
the results of the measurements are usually interpreted in 
terms of a local resistivity tensor. Due to high anisotropy 
of these materials the in-plane resistivity p xy is much 
smaller than the out-of-plane resistivity p z . Commonly 
the resistivity is assumed to be a local function of the 
current density, and to depend on the applied magnetic 
field and the temperature 2-4 . A recent experiment by 
Khaykovich et al. 1 does not fit into this scheme. In this 
experiment transport and magnetization measurements 
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FIG. 1. Resistance at various I a (left axis, log scale) 
and magnetically measured critical current (right axis, linear 
scale, open circles) vs T for the irradiated sample, H a = 200 
Oe (a) and H a = 500 Oe (b) (taken from Khaykovich et al. 1 ) 



in BSCCO crystals at elevated transport currents and 
perpendicular magnetic field are performed, using high 
quality BSCCO platelets with current leads attached to 
the top surface and an array of 2DEG Hall sensors to 
the bottom surface. At a first glance, the pictures that 
emerge from the transport and the magnetization mea- 
surements are mutually contradicting. Transport mea- 
surements reveal finite resistivity below the magnetic ir- 
reversibility line, in the superconducting state. This re- 
sistivity is non-monotonic with temperature, showing re- 
entrant behavior, and non-linear with current. As seen in 
the graphs of R vs. T, Fig. 1, at low transport currents 
R(T) is monotonic, dropping below experimental resolu- 
tion when temperature is reduced. At elevated currents, 
the resistance initially drops as T is lowered, but then 
goes up, the bump being steeper at lower currents. Also 
R{T) shows strong non-linearity, so that an increase of 
the current by 30% or less may result in enhancement of 
R by orders of magnitude. The source of this resistance 
is, presumably, vortex flow as a response to the electric 
current. 

In contrast, local magnetization measurements in the 
presence of transport current, shown in Fig. 2a, indicate 
that the vortices are pinned. These measurements can be 
well described in terms of the Bean model of the critical 
state 6,7 . The model states that below the irreversibil- 
ity line the local current density equals cither zero, or 
the critical current density, directed in such a way as 
to obtain the total transport current and the magnetiza- 
tion. The spatial distribution of the magnetic field is then 
given by the Biot-Savart law 8 . Since the current density 
nowhere exceeds the critical one, the Bean model predicts 
zero resistance. Within the Bean model finite resistivity 
can be expected only above the magnetically measured 
irreversibility line, which in Fig. 2a occurs above 1600 
Oe. Indeed, at low I a the measurements (carried below 
the irreversibility line) show practically zero resistance. 
However, at elevated currents, substantial resistance is 
measured concurrently with the hysteretic magnetization 
well below the irreversibility line, as seen in Fig. 2a. Fig- 
ure 2b shows the corresponding field profile B z {x), ob- 
tained by the array of Hall sensors at 400 Oe in presence 
of transport current on increasing and decreasing H a . A 
clear Bean profile is observed. Fitting this profile to the 
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FIG. 2. (a) Resistance (right axis) and hysteretic magneti- 
zation loop in the sample center (left axis) vs. H a at T=30 K 
and I a =25 mA. (b) Profile of magnetic induction across the 
sample at 400 Oe on increasing (•) and decreasing (o) fields 
(taken from Khaykovich et al. 1 ). 

theoretical field distribution in platelet sample results in 
total critical current of I c = 4.2 A, which is more than 
two orders of magnitude higher than the transport cur- 
rent of 25 mA. Figure 1 shows I C (T) determined from 
the Bean profiles together with the resistive data. The 
re-entrant resistance always occurs in the region where 
zero resistance is expected, since the transport current is 
much lower than the critical current. 

Thus, the main puzzling observations of Khaykovich et 
al. are the nonvanishing resistance below the irreversibil- 
ity line, which indicates flux flow, coexisting with mag- 
netization measurements which indicate that the vortices 
are pinned, the re-entrant behavior of the resistance with 
the temperature and its strongly nonlinear dependence 
on the current. 

Khaykovich et al. 1 suggest the following qualitative un- 
derstanding of the observation. BSCCO, being a strongly 
anisotropic type II high T c superconducting material, 
consists of superconducting CuC>2 layers, separated by in- 
sulating barriers. Each layer can carry current, resulting 
in total parallel current along the sample. Also, due to 
Josephson coupling between the layers, current can flow 
perpendicular to the layers. Because of large anisotropy 
a typical ratio of the perpendicular and parallel resis- 
tivities is ~ 10 4 in the normal state. In perpendicular 
magnetic field the flux penetrates the system in form 
of vortices, but, due to weak interlayer coupling, these 
are two dimensional "pancakes" , rather than three di- 
mensional filaments. Pancakes in the same layer repel 



one another, while those in different layers attract via 
Josephson and magnetic coupling 7 . In the experiment, 
the leads are attached to the top surface of the crystal. 
Hence the current distribution is non-homogeneous along 
the sample thickness, planes near the bottom of the crys- 
tal carrying much lower current than those at the top. 
As temperature decreases, pinning of vortices becomes 
more effective. Eventually the critical current density 
exceeds current density near the bottom. Then pancake 
vortices at the bottom stop moving, while pancakes at 
the top maintain their high velocity, since current den- 
sity there is much higher than the critical current den- 
sity. As a result, velocity gradient of pancake motion be- 
tween different layers is increased. This, in turn, leads to 
shear-induced phase slippage between the adjacent CuC>2 
planes, reducing the Josephson coupling and increasing 
the perpendicular resistance p z . The larger p z causes the 
current to flow in a thinner part of the sample, thus mak- 
ing the process self-enhancing. Since all of the transport 
current flows in a few layers near the top of the sample, 
finite resistance exists at currents much lower than the 
critical current expected from the Bean model. Magne- 
tization measurements, on the other hand, measure the 
magnetic response of all layers. When the vortices are 
pinned in most layers, this response is irreversible. 

In this work we take this qualitative explanation as 
a starting point and construct macroscopic and micro- 
scopic models to analyze the experiment. We start by 
presenting a macroscopic model in which the sample is 
assumed to be constructed of a resistive part, an inter- 
face and a dissipationlcss part. The perpendicular resis- 
tivity of the resistive part is assumed to depend on "vor- 
tex shear" . The parameters of this model are introduced 
phenomcnologically. We then examine the dependence of 
the sample's resistance on these parameters, and the con- 
clusions that may be drawn regarding the dependence of 
the resistance on the temperature and current. Following 
that we construct a microscopic model aimed at deriving 
an expression relating the conductivity in the direction 
perpendicular to the layers to the inter-layer variation 
of the current parallel to the layers. Finally we compare 
the conclusions of our model to the experimental findings. 
Although we find a general agreement, we also point out 
some remaining difficulties, associated mostly with the 
lack of quantitative information regarding several of the 
parameters of the model. 



II. THE MACROSCOPIC MODEL 

As we focus here on the consequences of inhomogene- 
ity in the current distribution in z direction, we use a 
one-dimensional model in which all quantities can vary 
only in this direction. Since scales of interest are much 
larger than the microscopic scale defined by the spac- 
ing between adjacent superconducting layers, we take a 
continuous limit in z direction. 
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FIG. 3. The macroscopic model. 

The model is described as follows. A current Jj„ is 
injected into a system of depth d from above. Part of 
the current then flows horizontally as j x and the rest 
- vertically down as J z (then, of course, returning ver- 
tically up at the other end of the system). Below the 
depinning temperature T<j, when there is non-zero criti- 
cal current density j c , the system can be divided into two 
parts, by the value of the in-plane current j x . The up- 
per part of the system carries current density larger than 
j c , so it has finite resistance, while at the lower part the 
current density is smaller than j c , and thus it has zero 
resistance. Accordingly, we consider the system as con- 
sisting of two phases: a resistive phase at the top, having 
parallel resistivity R x and perpendicular resistivity p z /2, 
and a dissipationless phase with zero resistivity. Note, 
that since the current first flows down and then up, the 
total perpendicular resistivity it experiences is p z . Fur- 
thermore, we assume that current crossing the interface 
between the two phases faces a resistance Ri n t/2. The 
position of the interface is determined by the condition 
jx = jc- This condition also fixes the current J ou t flowing 
through the dissipationless region: 



J outRint — jcRx- 



(1) 



At high temperatures j c is zero, and the system consists 
only of the dissipative phase. 

The basic equations governing the distribution of the 
current in the dissipative phase are the two Kirchoff equa- 
tions. The continuity equation is (note that in the geom- 
etry we consider J z and j x have different dimensions, 
since J z is a two dimensional current density, while j x is 
a three dimensional current density): 



d z J z + j x = 



and the equation giving the total voltage is: 

V= f J z (z') Pz (z')dz'+j x (z)R x . 
Jo 



(2) 



(3) 



As we show below in the microscopic analysis, the z-axis 
resistivity depends on the difference between j x in adja- 
cent layers d z j x , and this dependence may be approxi- 
mated by 



Pz = Po + s[p\ + (fd zJx ) 2 = Po + \j P \ + {fdlJ z )\ (4) 

while R x is assumed to be a constant parameter. The 
term fd z j x in the resistivity p z is a contribution of the 
" shear" between vortices in different layers to the out-of- 
plane resistance. It originates from the effect of a veloc- 
ity gradient between vortices in adjacent planes on the 
Joscphson coupling between the planes. 

Substituting Eq. (2) into Eq. (3) and differentiating 
with respect to z we obtain: 



Mpo + yffi + (fd 2 z J z f) - R x d 2 z J z = 0. (5) 
This equation can be solved only if the condition 

(6) 



J z < — d z J z 

PO 



is satisfied. Designating J = J z , J" = d z J z and solving 
for J", we obtain 



J" = 



1 



(Rx/jy p 

x (r xP o/j ± yWiA/) 2 + / 2 (pg-p?)) • (?) 

The condition (6) requires that plus sign be taken in 
Eq. (7) and that J < Jq = R x /.f- This means that 
as Ji n — ► Jo, both J" and J' diverge, so that the voltage 
V also diverges, and the system becomes insulating. In 
fact, as J in — > J , current gradients in the system become 
large, and then the quasi-particle channel for z-axis cur- 
rents needs to be taken into account, as analyzed below. 
When doing this, we find that J is actually not a cut- 
off value for the injected current, but rather a parameter 
that signifies the importance of shear effects. Thus, when 
J in becomes comparable with J , shear becomes strong, 
and the resistance is strongly non-linear with J in - 

Substituting the solution for J" into Eq. (4), the per- 
pendicular resistivity can be expressed in terms of J: 



Pz = Pa 



1 + v/r 2 +K 2 (l-r 2 ) 



(8) 



where we used reduced quantities n = J/ Jo and r = 
pi I po- It is plotted in Fig. 4. Again, this is valid for J 
not too close to J . 

It is possible to integrate Eq. (7). Some intuition to it 
may be obtained by noticing that Eq. (7) may be viewed 
as an equation of motion for a particle whose one dimen- 
sional coordinate is J, its "time" is z, and the potential 
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FIG. 4. Perpendicular resistivity p z for the model without 
quasiparticle channel, Eq. (8) (dashed line) and with quasi- 
particle channel, Eq. (12) (solid line). 



it moves in is 



U(J) = 



JoPo 



r 2 + k 2 (1 — r 2 ) 



+ ln|l - v/r 2 + «: 2 (l-r 2 )|), 



(9) 



This potential is plotted in Fig. 5. 

The analysis leading to Eq. (7) neglects inter-layer 
current flow by means of quasi-particlc tunneling. When 
the current gradient J" gets large, p z becomes large, and 
a large portion of the current flows perpendicularly in the 
form of quasiparticles. Hence in this high-gradient limit 
the perpendicular resistivity should be modeled by two 
resistors in parallel. Also, since in this regime the current 
gradients are large, a linearized expression for the Joseph- 
son channel resistivity, Eq. (4), can be used. Thus the 
Josephson channel carries a resistivity po + fd z j x , while 
the quasiparticle channel's resistivity is p qp . The total 
perpendicular resistivity is 



P z 1 {d z j x ) = (p + fd z j x ) 1 + P, 



qp 



(10) 



It is assumed, of course, that p qp po- Using this as- 
sumption and solving again for the current distribution, 
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FIG. 5. Potential U(J) (taken with minus sign) for the 
model without quasiparticle channel, Eq. (9) (dashed line) 
and with quasiparticle channel, Eq. (13) (solid line). 



we get a differential equation 

J " = ff ( K ~ 1 ± V( k - 1 ) 2+4k ^p) ' ( U ) 

where again k = J/Jq. In order to have J" > 0, we need 
to choose the plus sign. The perpendicular resistivity can 
be expressed in terms of J: 



Pqp' 



1 + y/( K - l) 2 + 4«r 2 
2k 



(12) 



where r 2 = Po/p qP - It is plotted in Fig. 4. 

Integrating Eq. (11), we obtain the corresponding "po- 
tential", plotted in Fig. 5: 



U(J) 



JoP, 



qp 



4/ 



(k - l) 2 + 4r 2 (l - r 2 )arcsh 



k - 1 + 2r 2 
2^2(1-7-2) 



-(«-! + 2r 2 )v/4r 2 (l - r 2 ) + (k - 1 + 2r 2 ) 2 



• (13) 



In both cases we may use U(J) together with the bound- 
ary conditions to determine the resistance of the system. 
The "velocity" of the particle is given by 



Z J Z = -y/2[C-U(J)]. 



(14) 



Here C is a constant determined by the boundary con- 
ditions, which require J z (z = 0) = Jj„, and either 
J z {d) = 0, in the case where the whole sample is resistive 
so that j x > j c \ or J z (d int ) = J out and d z J z (d int ) = j c , 
for the case where the lower part of the sample is dis- 
sipationless and is separated from the upper part by an 
interface at depth dint ■ The first case takes place at tem- 
peratures above Td, where j c = 0, while the second case 
- when the temperature is below Td- For the second case 
the derivative J'(z = 0) can be found: 



J'(0) = -^2{U{J out )-U{J m ))+ 3 l 



(15) 



Finally, after solving for the current profile J(z) we may 
calculate the resistance of the sample to be 



RxJ'jQ) 
j(0) • 



(16) 



We now summarize how the resistance depends on vari- 
ous parameters of the model. 



A. Dependence on the injected current J;„ 

For high temperatures, when the whole sample is dis- 
sipative, the resistance increases monotonously with Ji n , 
experiencing a sharp increase around Jq. This is because 
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larger currents produce larger current gradients, which, 
in turn, increase the vertical resistivity p z . 

The low temperature case, where there is an interface, 
is more complicated and depends on the value of the in- 
terface resistance Rint- To investigate the dependence of 
the resistance on the current, we need to differentiate the 
expression 



_ R x J'(0) = R x ^2{U(J out )-U{J in ))+f c 

Jin Jin 



R = - 

with respect to Ji n . We have: 



(17) 



dR 

J J in 



R x dJ'(O) J'(0) 
~j T7 h Hx ~ 

Jin O.Jin 

Rx 



J 2 

in 



J m J'(0) dJ, 



dU{J in ) , „ J'(0) 
— - • R, — . 



(18) 



Substituting dU{J)/dJ = J" and multiplying by a posi- 
tive quantity — R x J'(0), we obtain 



dR R 2 x J"(0) 

oc 



dJ i] 



J'ii 



R 2 (Jin)- 



(19) 



The above derivative is definitely positive at J^„ « Jo, 
since, as we saw above, the Cooper pair channel gets 
blocked, and the resistance of the system rises abruptly 
as J in approaches Jo- Hence it is only left to determine 
the dependence on J in for J in « Jo- For this case we 
can neglect the quasi-particle contribution and use the 
expression (7) for J"(0). After a minor manipulation we 
obtain 



dR 

dJ in. 



OC 



R xPo (l-r 2 ) 



\/r 2 + re; 



(1-r 2 ) 



R {Jin), 



(20) 



where re in = J in /J - 

The sign of this expression determines whether the re- 
sistance increases or decreases with the injected current 
Jin- It is easily verified that this expression is increasing 
with Ki, n (i-e. with Jj„). Hence it is enough to determine 
the sign at the smallest current at which the model is ap- 
plicable, J in — J out- if h is positive, then the resistance 
increases monotonically with the current, while if it is 
negative, the resistance first decreases and then starts to 
grow as the current becomes large enough, c.f. Fig. 6. 
Substituting re in = K out and R(J ou t) = Rint, we get: 



dR 



d J i' t 



Jin 



cx 



R xPa (l-r 2 



l-^r 2 + [jJ/R~tKl - r 2 ) 



-R 



int ■ 



(21) 



The result is a decreasing function of Ri n t ■ It is positive for small Ri nt (which should be larger than fj c in order to 
satisfy J out < Jo), negative for large Ri n t and vanishes at Ri n t = Rintfi, given by 



Rint,o = \j R.P0 + J 2 f 2 /2 + J (R xP0 + .f c f 2 /2) 2 - R 2 p 2 (l -^). 



(22) 



Thus, the dependence of the resistance on the current 
is controlled by the value of Rint, as is seen in Fig. 6. 
F° r fjc < Rint < Rint,o the resistance increases mono- 
tonically with the current. But if Ri n t > Rintfi, the 
resistance decreases for small currents J out < Ji n < Ji, 
where J\ is the solution of an equation 



RxpoQ 



r 2 ) 



1- V^ 2 + (JiM)) 2 (l-r 2 ) 



= R 2 (Ji) 



(23) 



The physical explanation for this behavior is that when 
the current is increased, the interface is pushed down- 
wards, increasing the thickness of the upper (dissipative) 
layer. If the interface is highly conducting (small Rint), 
most of the current is shunted through the lower (dis- 
sipationless) part of the system, so the increase in the 
upper layer thickness increases the resistance of the sys- 
tem. However, if the interface is almost insulating (large 
Rint), most of the current flows through the upper part, 
and by increasing its thickness the resistance of the sys- 
tem is decreased. Of course, at large enough currents the 
rapid increase of p z due to shear has dominant effect, so 
the resistance increases anyway. As we show below, the 



relevant case is large Rint, when a re-entrant behavior 
as a function of the temperature takes place. Hence, be- 
low Td an increase in the current influences the system 
in two opposite ways: it tends to decrease the resistance 
by moving the interface downwards; while through the 
effect of shear it tends to increase it. Also, we see that a 
strong increasing dependence of the resistance on the 
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FIG. 6. Sample resistance R{J in ) for R int < R int ,o (dashed 
line) and for 7?;„ t > Rintfi (solid line). 
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current appears only when Jj„ w Jo, this being true bot 
above and below Td- 

The difference between the results given here and the 
qualitative arguments of Rcf. [ 1] may be understood in 
the following qualitative way. Suppose that a current 
Ji n flows into the system and generates a current profile 
J(z) with an interface at z — d. When Jj„ is slightly 
increased one may expect the current gradient d z j x to 
increase, thus increasing p z , increasing anisotropy and 
pushing the interface upwards. The shear-induced in- 
crease in p z and the motion of the interface both tend 
then to increase the resistance. Our model yields a dif- 
ferent picture: as Jj„ is increased, the interface is shifted 
downward, thus reducing the resistance. The motion of 
the interface and the shear-induced increase of p z operate 
then in opposite directions. 



B. Dependence on the intralayer resistivity R x 

For the case when there is no interface in the system, 
increasing R x makes the current distribution more ho- 
mogeneous, so that j x (0) and j x (d) differ less. Put in 
another way, Jo = R x /f grows. Because of this, the ef- 
fects of inter-layer vortex shear become weaker, and the 
vertical resistivity p z decreases. Hence the total resis- 
tance R is influenced by two opposite effects: increase 
of R x directly increases R, this effect being dominant at 
small currents. On the other hand, through the decrease 
of p z it tends to decrease R, this effect becoming domi- 
nant at strong currents, when effect of shear is important. 
Hence the resistance grows with R x at small Jj n , while 
it decreases with R x as Ji n approaches Jo- 



C. Dependence on the critical current j c 

Next we discuss the dependence of the sample resis- 
tance on the critical current. We disregard a possible de- 
pendence of various parameters (like Ri n t , for instance) 
on j c and consider only a variation of R due to a shift in 
the position of the interface and the subsequent current 
redistribution. To determine the sign of the derivative 
J-p, we use the expression given by Eq. (17). First, it is 
easy to see that the sign of the derivative is independent 
of the value of injected current Jj n . Hence we can find it 
at Jm — Jout (he. when the interface is right at the top 
of the sample) . But when this condition is satisfied, the 
resistance of the system is constant and equal to Rint- 
Thus, 



djc 



dR dR_ 

djc dJ in 



dJ n 



Jin = Jo 



djc 



Then, using Eq. (1), we see that 

R x dR 



dR 

djc 



= 0. (24) 



(25) 
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FIG. 7. Sample resistance R(j c ) at Ji„ = 0.8 Jo 
for R int > W R x (po + Pi) (dashed line) and for 
Rint < \J Rx(po + pi) (solid line). 

Consequently, at small Jj„ the dependence of R on j c is 
opposite to its dependence on Jj„(at larger currents the 
dependence on j c remains of the same type, while the de- 
pendence on Ji n may change, as was shown above). This 
behavior is natural, as by increasing J in the interface is 
pushed downwards, while increasing j c it is pushed up- 
wards. Hence citing the previous results we obtain that 
R increases with j c for Ri„t > Rintfi, while it decreases 
with j c for Ri n t < Rint,o- Formulated in a different way, 
this means that R increases with j c when j c < j c o , where 
jco is given by 



JcO 



Ri nt l[l-R xPo {l-^)/RU 2 -r 2 



f 



l-r 2 



(26) 



In order that j c o be real and positive, Rint has to satisfy 
Rint > v Rx (po + Pi), which is physically plausible, as 
the interface should be insulating enough in comparison 
to the resistive phase in order that the rise in its verti- 
cal position would increase the sample resistance. If this 
condition is not satisfied, or if j c > j c0 , R decreases with 
j c - Note that when Ri n t — > oo also j c o — > oo, so that in 
this case R increases with j c for any relevant value of j c - 
The behavior of the resistance as a function of j c can be 
seen in Fig. 7. Also, from Eq. (17) we sec that ~ 
(since a square root is differentiated) . Hence at large cur- 
rents the dependence on j c becomes weaker. This is in 
contrast to the dependence on J in , which becomes very 
strong as Ji n — * Jo- All this, of course, is valid when j c 
is strong enough that there is an interface in the system. 



III. PERPENDICULAR RESISTIVITY DUE TO 
PARALLEL CURRENT GRADIENT 



The interplane transport properties of high-T c su- 
perconductors have been a subject of intense re- 
search over the past decade, both theoretical 9 ' 10 and 



G 



experimental 2 ' 4,11-13 . This transport, being of Joseph- 
son nature, is determined by the phase coherence be- 
tween the adjacent layers. For a superconductor in a 
perpendicular magnetic field, the pancake vortex struc- 
ture determines the above properties through the phase 
distribution. The vortex structure in high-T c supercon- 
ductors exhibits a rich variety of phenomena, including 
decoupling, melting, pinning, Bose glass formation etc., 
due to thermal fluctuations 14-17 , point defects 16-18 or 
columnar defects 17,19-22 . For the perpendicular resistiv- 
ity p z , the microscopic origin of the dissipation is less 
obvious than for the in-plane resistivity p x , where it is 
understood in terms of the Lorentz force, acting on the 
pancakes. Koshelev 10 proposed a microscopic mechanism 
for interplane dissipation, in which the pancake dynam- 
ics are shown to influence the interplane conduction, and 
calculated p z for the simplest case of non-interacting pan- 
cakes. Following Ref. 10 we analyze a simple microscopic 
model aimed at a derivation of a formula for a contri- 
bution to the resistivity in z direction (perpendicular to 
the layers) pl° of a superconducting slab due to a gradi- 
ent in the current in x direction (parallel to the layers). 
We first derive p\ D for a 3-dimensional sample assum- 
ing no interactions between the vortex pancakes. Then 
we show how the results are modified in presence of in- 
terlayer and intralayer correlations between the pancake 
positions. Finally, we transform the 3-D resistivity pa- 
rameters into a form appropriate for the 1-dimensional 
model used in the previous section. That is, we show 
how R x , p and / of the macroscopic model are derived 
from the resistivities of the 3-dimensional model. 



A. Noninteracting pancakes 

We assume a layered superconductor with noninter- 
acting pancake vortices in it. The vortices arc mobile, 
and their relative diffusive motion provides a mechanism 
for perpendicular resistance. In addition, each layer car- 
ries a different current, causing different drift velocities 
of vortices in adjacent layers. This increases the decay 
of phase correlations in time, thus enhancing the perpen- 
dicular resistance. 

We start from the Kubo formula for finite tempera- 
tures: 

al D = J dr dt (sin 5<j)(0, 0) sin d<f>(r, t)). (27) 

Here s is the interplane separation, jj - the Josephson 
current and S(f> - the gauge invariant phase difference be- 
tween neighboring layers. We neglect interplane correla- 
tions, so that averages like (expi<5</>) are assumed to be 
zero and then 

(sin(ty(0,0)sin(ty(r,i)) w (l/2)i?e(exp[iS(r,t)]>, 

where S(r, t) = 8<j>(r, t) - 50(0, 0). (28) 

Next we assume Gaussian randomness of S, so that 

(exp[iS(r,t)]) =exp[-<S(r,t) 2 )/2]. (29) 



Thus we need to calculate the mean square of S. We 
write 

S(t, t)=^<t> v {r- Ri,i(t)) - <Mr - R 2 ,i(t)) 

i 

-0„(-Ri, j (O)) + ^(-R 2 , i (O)), (30) 

where <^(r) is the phase distribution of a single vortex. 
Expanding, we write 

S(r, t) = [r - AR M (t)]V<M-R M ) 

i 

-[r-AR 2ii (t)]Vi(-R 2 ,i), (31) 
where AR(i) = R(i) - R(0) and 

i x r 



V<£„(r) = 



(32) 



Now we assume that the pancakes in the layers are ran- 
domly placed, so that 

(V0„(-Ri)V0„(-R2)) =0 and 
(V0„(-Ri,i)V^(-Rij)) =MV0„(-Ri,i) 2 >. (33) 

Then the square of a sum breaks into a sum of squares, 
so that 

(S(r,t) 2 ) = £ <([r-AR M (i)]VM-Ri,i)) 2 > 

i 

+ (([r - AR 2 ,,(t)]Vi(-R 2! ,)) 2 ). (34) 

Now we can write for each layer 

AR(t) w vt + SR(t), (35) 

where v is the drift velocity of vortices due to the cur- 
rent, and 5R(t) is the diffusion term. It gives the main 
contribution at zero current gradient, and we will copy 
it from the Koshelev's article. Using the expression (32) 
we write 



(S(r,i) 2 )=£(r- Vl t) 2 



+ (r-v 2 *) 2 < 



R 



Ri 



R 2,i 



) + (S diff (t) 2 }. (36) 



Now we calculate the averages: 



dR 



■ In 



Rr, 



Rmin 



(37) 



where n is the density of the vortices and R m i n and R m ax 
- the lower and upper cutoff radii. Substituting this, we 
obtain 
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<S(r,i) 2 ) = [(r - v lt f + (r - v 2 i) 2 ]™ In §^ + <S»„(t)> 

i? 



= [2(r - Vt) 2 + (Avt) 2 /2]7rnln + (S^ (*)>, (38) 



L min 

where 

V = (vi+v 2 )/2 and Av = vi-v 2 . (39) 
Substituting this result back into Eq. (27) and using Koshelev's result for Sdiff, we obtain: 

a 3 z D (Av) =*£\ drdtexp f-[(r - Vi) 2 + (Avt) 2 /4]™ In - 2™^ ln(i? 2 / R 2 mm ) 



Sjj 1 



2T nln(i?j/a ) Av^im ln(Rj/a a ) 



F (4D A /7rnln(/2j/ao)/Au) , (40) 



where D is the diffusion constant of pancake motion inside the layers. We used the Josephson radius Rj for the upper 
cutoff radius and the average intervortex spacing ao - for the lower cutoff. The function F(y) is defined by 



F(y)= dxe~ x2 - 2xy = e y2 / dxe~ x ' = ^—e^ [1 - Erf(y)] . (41) 
This function can be easily approximated for small and large values of its argument: 

™-{^-? < 42 » 

Using this and expressing the vortex velocity difference in terms of parallel current gradient modulus, 

Av^^ S 2 (<P /c)d z f x D , (43) 

where the average pancake mobility \i is connected with the diffusion constant D by the Einstein relation D = /j,T, 
we obtain for the perpendicular resistance 



o 3D (d .am = 1 2^ 5 ($ /c)a z j^[nln(i? J /q )] 3 / 2 exp(-7rnln(i 1 ', / /a )[4r/ S 2 (cI> / c )9 2 j^] 2 ) 
3 j 1 - Eii[4T^nln(Rj/a )/s^ /c)d z jl D } 



Expanding this, we obtain for small current gradients: 



pl D (d z .f x D ) = ^ (8nD[nln(Rj/a )} 2 + ± n ln(Rj/a )D[ S -^-d z j 3 x D } 2 j , (45) 
i.e. a parabolic dependence on d z j 3D . On the other hand, for large current gradients, 

P 3 z D (d z f x D ) = ™ ({n\n(Rj/a )}V 2 S 2 (<t>o/c)d z .f x D + 8[nln(Rj/a )} 2 T) (46) 



i.e., a linear dependence on d z j 3 



3D 



As Eq. (44) is not convenient for analytical work, ^3£> _ 16-Pr ^ m(_Rj/a )] 2 

we will use an approximation of the form p 3D {d z j 3D ) = sj 2 

PI D + V(PI U ) 2 + {P D d z f x D ) 2 which gives a correct value 3D = m EL\ nUR /a )f 

at zero current gradient and the asymptotic behavior at Pi \ ) g j2 i ( J i °n 

large current gradients. It also approximates quite well 2Ds $ 

the behavior of p 3D (d z j 3D ) in the intermediate range of f 3D = — [n\n(Rj / ao)} 3 ' 2 ■ (47) 

current gradients. Comparing the coefficients, we obtain - 7 J c 
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B. Correlations between pancake positions 

Here we demonstrate how the results obtained above 
are modified in presence of inter- and intralayer correla- 
tions between pancake positions. 

We first consider the effect of interlayer correlations. 
The presence of such correlations can be crudely de- 
scribed by regarding pancakes in different layers as tied 
together into vertical line segments of length L z , which 
move as a whole. These segments should be used instead 
of independent pancakes of previous subsection. The 
phase differences 6(j>(r,t) and corresponding Josephson 
currents are created only at the ends of these segments 
(more exactly, between layers, where one segment ends 
and another one starts), while the middle parts of the 
segments do not contribute to 6<j)(r,i). This means, that 
the effective concentration of vortices is reduced by a fac- 
tor L z /s. Next, since each line segment has an increased 
"mass" , the mobility fj, and the diffusion constant D are 
now reduced by another factor L z /s. Finally, the ver- 
tical separation between the segments is L z instead of s 
for free pancakes. This means that the velocity difference 
between the segments due to current gradient is increased 
by L z /s. To take into account this and the reduction in 
the mobility in Eq.(43), the flux quantum $0 should be 
multiplied by (L z /s) 2 . This specifies, how the resistiv- 
ity parameters are modified in the presence of interlayer 
correlations. 

Next we turn to consider the intralayer correlations. 
Roughly speaking, these correlations cause pancakes in 
each layer to aggregate in clusters of size L xy , so that 
there are (L xy /a ) 2 pancakes in a cluster. Pancakes 
inside each cluster are ordered, while different clusters 
move independently (actually, there is a hard-core repul- 
sion between them). Since vortices in the same cluster 
are not independent, Eq.(34) for the phase correlation 
square now reads as 

<S(r, tf) = £<(]>>- AR M (t)]V^(-R M )) ) 

+([^[r-AR 2ii (t)]Vi(-R 2ii )] ), (48) 

\iea / 

where a is an index of a cluster, while i - of an indi- 
vidual pancake. For clusters which are far enough away, 
the differences in the location of individual pancakes in- 
side the cluster can be neglected. Then each such clus- 
ter gives a contribution to (S(r,t) 2 ), which is {L xy /a ) 4: 
times larger than a contribution of an individual pancake. 
On the other hand, the concentration of the clusters is 
n(ao/ L xy ) 2 . To take both effects into account, we should 
multiply n by (L xy /ao) 2 in the final result. Also, the 
diffusion constant (and the mobility) of each cluster is 
reduced by a factor (L xy /a n ) 2 7 while the flux quantum 
$0 should be multiplied by the same factor. Substituting 



all these prescriptions into Eq.(47), we obtain the resis- 
tivity parameters in the presence of correlations between 
pancake positions: 

3D WDT (L xy /a ) 2 2 

P l» = (87r _ 16) ^i^y/^ [nHRj/ao) f 

S J.J \ L zl s ) 

f D ^^i^^ [nHRj/aQ)] m. (49 ) 

Here we neglected all changes in the argument of the log- 
arithms. 



C. Transformation of the parameters into ID form 

Now we transform these quantities into a form appro- 
priate for the ID macroscopic model. For this, we first 
define the corresponding fields and currents from their 
3D counterparts (assuming that everything is uniform in 
y direction): 

j x (z)=j 3 x D (x = L x /2,z)L y 
J z (z)= / dxj 3D {x,z)L y 



V(z)= C' dxE 3 x D (x,z) 
Jo 

E z (z) = E 3 z D (x = 0,z), 



(50) 



where L x and L y are sizes of the sample. Then, for large 
current gradients, we use the Ohm's law for the 3D sam- 
ple and average over x: 



/ * dxE 3D (x,z) = f * dxpl D 3 3D (x,z) 
Jo Jo 

f L x /2 rL x /2 

j o dx E 3D (x,z) = J dx [p 3D j 3 z D (x, z) 

+f 3D d z j 3 x D (x,z)j 3D (x,zj). (51) 

Defining now the reduced quantities as ratios between 
the 3D and ID ones, so that 



3l e \x,z), 



3x{z) 



Aredi \_ 3z D ( x i z ) 

Jz [X,Z) - J z (z) 

■red/^ \ E x (x,z) 



ET{x,z) 
E? d {x,z) 



V{x,z) 
E 3D (x,z) 



E z {z) ' 

we obtain from the previous equations: 



(52) 
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V(z) dxE™ d {x,z)=j x {z)pl D f Lx/2 dxj^ d (x,z) 

Jo Jo 

E z {z) ( L ' J2 dxE^ d (x,z) = J z (z)pl D f L ' 
Jo Jo 



+ Jz{z)f 



3D 



dxf z ed {x,z) 
L ' J2 dxf z ed (x,z)d z j x (z)f x ed (x,z). 



(53) 



Then, in order to obtain the equations of the macroscopic model, we make two assumptions: first, we neglect the 
derivative d z j x ed (x, z); second, we assume that the reduced quantities are not affected by shear effects, so we calculate 
them from a linear model with / = 0. The parameters of the macroscopic ID model are then given by 



R, 



dxjr d (x, Z ) 



= pI d I 

Jo 

Jo dXJ * ed{X ' Z) I J, 

-L x /2 



W 2 



dxE r x ed (x,z) 



po/2 = pl D j LX ' 2 dxf z ed (x, z) / [ LJ2 dxE™ d (x, z) 



//2 = / 



/ 

Jo 



dxf z ed (x,z)f x ed (x,z) 



' rL x /2 

Jo ' 



dxE r z ed ( Xl z), 



(54) 



where po and / are divided by 2, since, as we explained in 
the beginning of the previous section, the perpendicular 
resistivity of the macroscopic model is taken to be p z /2. 

To find the reduced quantities, we need to find the cur- 
rent distribution in a sample with constant resistivities 
p x D and pl D . This amounts to solving the Laplace equa- 
tion with the boundary conditions j x (x = 0, z) = j x (x = 
L x ,z) — j z (x,z — oo) = 0, and j z (x,z = 0) = 0, ex- 
cept two narrow regions near x = and x = L Xl where 
j z (x, z = 0) is, respectively, positive and negative. This 
describes contacts, attached to the top of the sample, 
where the current flows into and out of the system. For 
simplicity we assumed here that the system is infinitely 



thick in z direction. Choosing an appropriate form for 
j z (x, z = 0), we obtain: 



sinh k (w + az) cos k x 
sinh 2 k (w + az) + sin 2 k^x 



n sin kox cosh ko{w + az) 

Jx = «" 



(55) 



sinh 2 ko{w + az) + sin 2 kox ' 

where w is the width of the contacts (w <C L x ), a = 
V Pz° / Px D i an d = k/L x . Using this, we calculate the 
integrals of the reduced quantities and substitute them 
into Eq. (54), thus obtaining 



Rx = P X D — 7— cosh k (w + az) log coth ko (w + az)/2 
irL,. 



Po/2 - pl D 
f/2 = f D 



7T 



f 



L x L y sinhfc (w + az) arctanl/ sinh ko(w + az) 

7T 1 

2L X L^ sinh 2 ko(w + az) arctan 2 1/ sinhfco(w + az) 



(56) 



Here a finite z should be taken, so that k^az oc 1. Then 
the hyperbolic functions give factors of order 1, and the 
ID parameters are given by 



This establishes a correspondence between 3-dimensional 
resistivity parameters and the 1-dimcnsional ones, which 
were used in the macroscopic model. 



R x =pl°^ (57) 



Po/2 = pi 



nLy 

3D 7r 



IV. DISCUSSION AND CONCLUSIONS 



L x L y It is not easy to compare directly predictions of our 



/• /o — f 3D 77 model with the experimental results, since we do not 

2L x L y ' know temperature dependence of various parameters of 
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the model. Hence we make only qualitative statements 
based on robust features of the model. 

First, the model predicts that the resistance grows with 
the current (at least for not too small currents), and this 
current non-linearity becomes very strong as Jj„ — ► Jo- 
This is consistent with the experimental result. Using 
the results of the microscopic calculation Sec. Ill, we 
found that without correlations between the pancakes Jo 
is much larger than the relevant Ji n . However, in pres- 
ence of correlations its value is suppressed by a factor of 
(L xy /a )(L z /s) 5 / 2 , thus making its value much closer to 
Ji n . If the ratios L xy /ao, L z /s are assumed to be 10-15, 
Jo becomes comparable with the experimentally relevant 
currents. This provides an explanation to the experimen- 
tal fact that the current non-linearity becomes strong be- 
low the depinning transition temperature Tj, where cor- 
relations between the pancakes start to build up. Next, 
the model explains the feature of re-entrance, that is, 
the experimental observation that below the depinning 
transition the resistance increases as the temperature is 
decreased. According to the model, if the interface resis- 
tance R int is large enough, the resistance of the system 
grows with j c , which naturally starts to grow as the tem- 
perature is decreased below X^. Moreover, the model 
predicts that this rise in the resistance should be more 
pronounced for smaller currents, as indeed observed. 

Some ingredients are missing from our model. First, 
the model approximates R x to be independent of the 
intra-layer current. This approximation is presumably 
good above the depinning temperature, but becomes 
poor below that temperature, where intra-layer current 
induces vortex depinning. Second, a missing ingredient 
in our work is a microscopic derivation of the interface 
resistance Rint , separating between the resistive and non- 
rcsistive parts of the sample. The microscopic origin we 
have in mind is that in the region between the two phases 
the pancake mobility is very sensitive to parallel current 
variation. Then a small current gradient is enough to cre- 
ate a large pancake velocity gradient, which would cause 
a large perpendicular resistance in that region. Our at- 
tempts to provide a microscopic derivation of Rint and 
its temperature dependence led us to results that heavily 
depend on various microscopic parameters whose values 
and temperature dependences are not known. We were 
therefore led to leave R in t as a phenomenological param- 
eter. 

Altogether, then, our work is able to explain the qual- 
itative features of the non-linear transport observed in 
Rcf. 1 and unravel a unique feature of transport in super- 



conducting BSCCO samples in perpendicular magnetic 
field. 
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